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i We numerically study a nondisordered lattice spin system with a first order liquid-crystal tran- 

. sition, as a model for supercooled liquids and glasses. Below the melting temperature the system 

can be kept in the metastable liquid phase, and it displays a dynamic phenomenology analogous to 
fragile supercooled liquids, with stretched exponential relaxation, power law increase of the relax- 
ation time and high fragility index. At an effective spinodal temperature T sp the relaxation time 
exceeds the crystal nucleation time, and the supercooled liquid loses stability. Below T sp liquid 
properties cannot be extrapolated, in line with Kauzmann's scenario of a 'lower metastability limit' 
of supercooled liquids as a solution of Kauzmann's paradox. The off-equilibrium dynamics below 
T sp corresponds to fast nucleation of small, but stable, crystal droplets, followed by extremely slow 
growth, due to the presence of pinning energy barriers. In the early time region, which is longer the 
O ' lower the temperature, this crystal-growth phase is indistinguishable from an off-equilibrium glass, 

both from a structural and a dynamical point of view: crystal growth has not advanced enough to 
be structurally detectable, and a violation of the fluctuation-dissipation theorem (FDT) typical of 
structural glasses is observed. On the other hand, for longer times crystallization reaches a thresh- 
old beyond which crystal domains are easily identified, and FDT violation becomes compatible with 
^ ■ ordinary domain growth. 
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I. INTRODUCTION 



s 

When a liquid is cooled fast enough below its melting point, crystallization is avoided and the system enters the 
supercooled phase. Relaxation time increases rapidly in this temperature regime, and when it becomes comparable 
with the largest experimentally accessible time the system falls out of equilibrium, remaining stuck in a disordered 
phase called structural glass [1]. 

The configurational disorder of a structural glass is not caused by the presence of intrinsic disorder in the Hamil- 
tonian, since this is just the sum of deterministic interaction elements. This notwithstanding, a lot of attention has 
been devoted to the phenomenological analogies between structural glasses and mean-field spin-glasses, which are 
systems directly containing quenched disorder in the Hamiltonian in the form of random couplings among the spins 
[2]. In particular, models of spin-glasses with p-body interactions and p > 3, seem indeed to have many features 
in common with structural glasses [3], with the obvious (but crucial) difference that in mean-field systems there is 
a purely dynamical transition which is absent in real finite-dimensional glasses. Despite the close phenomenology, 
however, it is natural to ask why we use models with quenched disorder to understand the behaviour of physical 
systems with deterministic Hamiltonians. 

Disorder, one of the characteristics of the glassy state, cannot be expected to arise spontaneously from a deterministic 
mean-field Hamiltonian, because in this case the mechanically stable configurations are ordered: the local field is the 
I ' same for all spins, and thus all tend to point in the same direction. Thus if we wish to stick to mean-field and the 
associated analytical advantages, we must resort to quenched disorder to produce glassy behaviour. The argument 
does not apply to finite range systems, however, and we can ask whether finite dimensional spin systems with no 
quenched disorder and p > 3 display a glassy phenomenology similar to structural glasses (we of course know that 
p = 2, i.e. the Ising model, is not glassy under many respects). This is a relevant question, first because it is important 
to understand what the key ingredients responsible for glassy dynamics are, and second because lattice spin systems 
are normally easier to study than off-lattice liquid models. 

The answer to this question seems to be affirmative. A number of spin models without quenched disorder now 
exist that reproduce a phenomenology similar to supercooled liquids and glasses. Generally, one can say that in all 
these systems the non-trivial dynamic behaviour is related to the presence of frustration, due to complex even if non- 
disordered interaction terms in the Hamiltonian. In the SS model of ref . [4], this is achieved by mixing nearest- and 
next-nearest-neighbor interactions, which gives an interesting phenomenology, although not precisely glassy. In the 
plaquette (PQ) model of [5,6,14], on the other hand, glassiness is achieved by a p-body interaction term with p = 4, 
in close analogy with p-spin models of mean-field spin-glasses. In the lattice glass model of [9] frustration arises from 
a constraint on the numbers of neighboring particles on the lattice. The PQ model has been heavily studied in recent 
years, and the results seem to suggest that its phenomenology has at least some points in common with structural 
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glasses. In particular, a sufficiently rapid cooling of the PQ model forms an off-equilibrium disordered phase which 
closely resembles real glasses. Therefore, it would seem that p-body interactions in finite dimensional spin models 
without quenched disorder are sufficient to reproduce a phenomenology similar to supercooled liquids and glasses. 

The aim of this paper is to study a model in the same class as those just described, namely the deterministic version 
of the CTLS model introduced in [7] (a concise description of some of our results has recently appeared in [8] ) . As 
discussed below, this model has various advantages with respect to its predecessors and will enable us to check more 
deeply to what extent lattice spin systems can mimic structural glasses. Moreover, important questions on the low 
temperature behaviour of glasses, such as the stability of the supercooled phase and the role of the crystal, can be 
addressed, whereas in previous models such an analysis was very limited. Let us now explain why. 

Liquids in nature have a crystalline ground state, usually with low degeneracy, and display a first order thermo- 
dynamic transition at the melting temperature. In the PQ model in three dimensions a first order transition takes 
place, but the 'crystalline' ground state has an unusual property for liquids, in that it is strongly degenerate: for 
linear size L, the PQ model has 2 3L different ground states with the same energy. The problem is not the degeneracy 
in itself (note that the ground state entropy density is zero) , but rather the consequences that this degeneracy has on 
the possibility to measure the degree of order in the system. More precisely, in the PQ model, given a ground state 
configuration (for example the ferromagnetic one), we can obtain another ground state by flipping any plane of spins. 
Clearly, by iterating this procedure we can produce ground state configurations which look highly disordered, despite 
all having the ground state energy. The problem then is how to distinguish a truly disordered, or glassy, configuration, 
from a configuration made up of many droplets of different ground states. The energy of both configurations will be 
high, in the first case because it is essentially a liquid, in the second case because of the interfaces among different 
droplets. In other words, energy is not a sufficiently precise measure of order in the PQ model, and due to the strong 
ground state degeneracy no configurational means can be devised to quantify order and domain growth. 

The system we introduce here does not suffer from this problem: it has a doubly degenerate ground state, such 
that it is possible to measure the amount of crystalline order and to monitor domain growth. A likely objection is 
that glassy physics does not have anything to do with the crystal, and that the possibility to measure the formation 
and growth of crystal domains is far from essential to study the physics of glasses. We reply to this objection with 
two observations. First, that simple glass formers usually do have a crystalline phase, whose possible connection with 
the properties of the glass cannot be completely ruled out in principle. In fact, our model explicitly shows such a 
connection. Secondly, if on the other hand spin models with a crystal phase fail somehow to reproduce some basic 
property of glasses, then we should understand why, and make sure that all spin models are free from this problem. 
Indeed, the precise nature of what is normally called the 'glassy' phase in systems like PQ [5] and in lattice glass 
models [9] is the main focus of our work. 

Our model shows, in this respect, an interesting phenomenology. As we will show, it has a first order liquid-solid 
transition at a melting temperature T m . Crystallization can be avoided by fast cooling, and a supercooled liquid 
phase, mctastable with respect to the crystal, is found. "Equilibrium" measurements can be performed as long as the 
equilibration time of the supercooled liquid is much shorter than the nucleation time of the crystal. In this phase the 
behaviour of our system is very similar to real supercooled liquids, and in particular to fragile systems. However, a 
temperature T sp exists, where the crystal nucleation time becomes of the same order as the liquid equilibration time: 
at this point the supercooled liquid looses stability and it no longer exists below T sp . We call T sp the effective spinodal 
temperature. It is important to stress the following fact: we are able to detect the loss of stability of the liquid only 
because the equilibration time at T sp , t sp , is shorter than our longest experimental time. We can state this point in a 
different way. A system falls out of equilibrium when the relaxation time becomes too large for the (real/numerical) 
experiments, and this is the operational definition of the glass transition temperature T g . Therefore, the reason why 
we can observe the loss of stability of the supercooled liquid in our model is that T sp > T g . The Kauzmann paradox 
(i.e. the fact that the extrapolated liquid entropy becomes smaller than the crystal entropy at a finite temperature) is 
avoided by this loss of stability: no equilibrium supercooled liquid exists below T sp and the entropy extrapolation is 
meaningless. In fact, this is the scenario proposed by Kauzmann himself [10] as solution of the paradox, which is thus 
exactly reproduced in this system. Below T sp the system enters a completely different phase, where crystal nucleation 
and crystalline domain growth are the main dynamic processes. However, by quenching to lower temperatures the 
growth of crystalline order becomes so slow that it is impossible to distinguish this phase from a typical off-equilibrium 
realistic glass. Moreover, we shall show that dynamical measurements of the fluctuation-dissipation (FD) ratio are 
compatible with the results obtained in ordinary structural glasses below T g . 

It is our impression that the scenario described above is also shared by the PQ, and possibly other spin models. In 
particular, that what is normally called the glassy phase for the PQ model is an off-equilibrium regime of very slow 
crystal growth. However, as already mentioned, without a direct measure of crystalline order it is hard to assess the 
amount of order in the PQ model, so the question whether there is in fact a glassy phase qualitatively different from 
the crystal growth phase below T sp is open. This is certainly not the case in the present model, though. 

Throughout the paper we will use a notation common to supercooled liquids and glasses. We call T m the melting 



2 



temperature. In analogy with mode coupling theory [11], T c will be the temperature where a power law fit locates 
a divergence of the relaxation time, while To will be the temperature where a Vogel-Fulchner-Tamman fit puts the 
same divergence. The so-called Kauzmann temperature, i.e. the temperature where the extrapolated liquid entropy 
becomes equal to the crystal entropy, is called T s . The effective spinodal temperature, i.e. the temperature below 
which it becomes impossible to observe the metastable liquid is called T sp . This temperature is called T g in [5], but 
we reject this choice because, as already mentioned, in the liquid literature the glass transition is defined rather as 
the point where the equilibration time is longer than the experimental time. 

The paper is structured as follows. In Section II we introduce our model, and describe the properties of the 
crystalline ground state. In Section III we show, by comparing the liquid and crystal free energies, that the model 
has a first order phase transition at a melting temperature T m . The dynamics of the supercooled liquid below T m 
and the metastability limit T sp arc considered in Section IV, while aging dynamics and its relation to domain growth 
are studied in Section V. Finally, we draw our conclusions in Section VI. 



II. THE MODEL 

The model is described by the Hamiltonian 

N 

h = J2m i + s *)> (!) 

i=l 

where the variables Sj = ±1 are Ising spins belonging to a two-dimensional square lattice of side L, and the plaqucttc 
variable is the product of the four nearest neighbors of spin i: 

h = ^sfsfsf , (2) 

(W is for west, S for south and so on). The model is studied by Metropolis Monte Carlo simulations with singlc-spin- 
fiip dynamics for lattices of size L = 100 and L = 500. 

A disordered version of this model was first introduced in ref. [7], as a way to model two-level systems in a lattice, 
interacting with their nearest neighbors, whence its name CTLS (Coupled Two Level Systems). In fact, one can take 
the spin Si to be the variable describing the state of a two- level system sitting at site i. Then state s» = — 1 has 
always zero energy, while the other state can be at 1 or -1 depending on the factor /j. Thus when one of the neighbors 
switches, the relative height of the two states at i is reversed. In [7], however, quenched random coupling constants 
where considered. Here we shall see that the multi-body nature of the interaction is sufficient to give a non-trivial 
dynamic and static behaviour. 



• 


o 





o 


o 


• 


o o o • o 


o 


o 


o 


• 








o 


o 


o 


o 


o 


o 


• 








o 


o • o o 


o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


• 


o 


o 


• o o o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


• 


o 


o 


o 


o 


• 





o 


o o • o 





o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


• o o • 


o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


o 


• 


o 


o 


o 





• 


o o o 





• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


• 


o 


• o o 


o 


o 


• 


o 


o 


o 


o 


• 


• 


o 


o 


o 


o 


• 


o 


o 





o • o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


o 


• 


o 


o 


o 


• o o 


• 


o 


o 


o 


o 


• 


o 


o 





o 


• 


o 


o 


o 


o 


• 





o o o 


• 


o 


o 


o 


o 


• 


o 


o 


• 


o 


o 


o 


o 


• • o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


o 


o 


o 


• 


o 


o 








• o o 


o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


• 


o 


• o o 


o 


• 


o 


o 


o 


o 


• 








• 


o 


o 


o 


o 


• 


o 





o o • 








o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


• o o 


• 


o 


o 


o 


o 


• 


o 


o 


o 





o 


o 


• 


o 


o 


o 





• 


o o o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


• 





• o 


o 


o 


o 


• 


o 


o 


o 


o 


• 


• 


o 


o 


o 


o 


• 


o 


o 


o 


o • o 


o 


o 


o 


• 


o 


o 


o 


o 





o 


• 


o 





o 


o • o o o 





• 





o 








• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


• 





o o o 


• 


o 


o 


o 


o 


• 


o 


o 


• 








o 


o 


• 


o o o o • 














• 











o 


o 


o 


o 


o 


• 


o 


o 


o 


o • o o 


o 


o 


• 


o 


o 


o 


o 


• 








o 


• 








o o • o o 








• 











o 


• 


o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o o • 


o 


o 





o 


• 





o 


o 





• 





o 


o 


o 


• o o o o 


• 





o 


o 





• 





o 


o 


o 


o 


o 


• 


o 


o 


o 


o 


• 


o o o 


o 


• 


o 


o 


o 


o 


• 


o 





o 


o 


o 


• 


o o o o • o o 





o 


• 








o 


o 


• 


• 


o 





o 


o 


• 


o 


o 


o 


o • o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


o 


• 


o 





O O • o o o 


• 


o 


o 


o 


o 


• 


o 


o 





o 


• 


o 


o 


o 


o 


• 


o o o o • 


o 


o 


o 


o 


• 


o 


o 


• 


o 


o 


o 


o 


• o • o 


o 


o 


o 


• 


o 


o 


o 


o 





o 


o 


o 


• 


o 


o 





o • o o o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


• 


o 


• o o 


o 


• 


o 


o 


o 


o 


• 


o 





• 


o 


o 


o 


o 


• 


o 





o o • o 





o 


o 


• 


o 


o 


o 


o 


• 


o 


o 


o 


• o • 


o 


o 


o 


o 


• 


o 


o 


o 





o 


o 


• 


o 


o 


o 


o 


• 


• 


o 


o 


o 


o 


• 


o 


o 


o 


o 


o 


• 


o 


• o o 


o 


o 


• 


o 


o 


o 


o 


• 


• 


o 


o 


o 


o 


• 


o 


o 





o • o 








o 


• 


o 


o 


o 


o 


o 


o 


• 


o 


o 


o 


• o o o 


• 


o 


o 


o 


o 


• 


o 


o 





o 


• 


o 


o 


o 


o 


• 





o o o • 





o 


o 


o 


• 


o 


o 



FIG. 1. Dextrocrystal (left) and levocrystal (right): empty circles correspond to positive spins, full circles to negative spins. 
In the left figure the square marks the five-element unit cell. In the right figure the small square shows the elementary unit 
used in the definition of the crystal mass (see sect. III). 



Frustration in the CTLS comes from the fact that it is not possible to satisfy simultaneously the 4-spin and the 
5-spin interaction. After careful inspection of the Hamiltonian it is found that the ground state of the CTLS is 
obtained by covering the lattice with the following non-overlapping elements: 

« = -l, *? = sf = *? = *? = +!■ (3) 
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As Fig. 1 shows, the covering can be done in two ways, giving rise to two different ground states which we call 
dextrocrystal and levocrystal. These two crystals are connected by the symmetry x — > —x, y — > y, or x — > x, y — > — j/, 
while they are both invariant under a; — > —a;, y — » —y. Therefore, the CTLS has a crystalline ground state, with 
energy density Eqs = —1-6. In Fig. 1 it can be seen that the crystal can be obtained by periodic repetition of a 
5x5 unit cell. Apart from the two distinct dextro- and levo- forms, additional crystals can be formed by translating 
the unit cell up to 5 lattice spacings in each direction. There are then different ways to cover the lattice, which are 
not locally compatible with each other. Upon a quench, the system will greedily optimize the energy locally in an 
uncoordinated way, creating many different competing domains with boundaries that turn out to be pinned, thus 
giving rise to slow dynamics. 

As already anticipated, the low degeneracy of the crystalline ground state is one of the main advantages of the 
CTLS compared to other plaquette models [5]. Indeed, as we shall see, in the CTLS it is possible to quantify the 
amount of crystalline order in the system by purely configurational means, irrespective of the value of the energy. 
This is particularly useful when we are studying the growth of crystalline domains, because in this case the energy 
may be very large, due to interfaces, even though there is a substantial amount of crystalline order. Moreover, the 
low degeneracy of the crystal makes this model more similar to real liquids. 



III. THERMODYNAMICS 



In the present section we will show that the CTLS has a first order phase transition at a melting temperature T m . 
To prove this fact we start from the numerical observation that in a certain interval of temperature both the crystalline 
phase obtained by heating the ground state and the liquid phase obtained by cooling a high temperature configuration 
are stable up to times of order 10 6 Montecarlo steps (MCS), for L = 100. In Fig. 2 we plot the equilibrium internal 
energy as a function of the temperature for the crystal and the liquid phases. 




FIG. 2. Liquid and crystal energy per spin as a function of temperature. The continuous lines correspond to the two fits 
reported in the text. The horizontal dotted line marks the energy of the crystal at zero temperature Ecs = —1.6. The vertical 
dotted lines marks T m — 1.30. L = 100. 



In order to find T m we need to compute the free energy in the two phases by integration of the internal energy To 
this aim, we fit the liquid energy (LQ) as 

E L q(T) = -a tanh(>/T), a = 1.84, 6=1.10, (4) 

while for the crystal (CR) we find 

Ec R (T)=E G s + cT n , c = 5.5xl0" 5 , n = 13. (5) 

By using the thermodynamic formula, 

0F(0) = 0oF(0o) + [ d(3'E((3'), (i = 1/T, (6) 

J A, 
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we obtain the free energy of the liquid and crystal phases: 

Flq(T) = -Tlog2 - £riogcosh(6/T), (7) 

F CR (T) = E GS - — t ——T n , (8) 
n — 1 

where we have taken /3 = for the liquid, and /3o = oo for the crystal. The melting temperature is fixed by requiring 
F L Q{T m ) = F CR (T m ), which gives 

T m = 1.30. (9) 

We have checked that our energy data are independent of system size for 100 < L < 500. 

In order to check that a first order liquid-crystal transition actually takes place at T m we proceed as in [5]: we 
consider an initial condition where half of the system is in the crystal phase and half in a random phase, and we 
run this configuration at T m + 0.01 and at T m — 0.01. If our determination of T m is correct, after a transient time 
the system should relax either to the liquid or crystalline phase, since the stable phase will expand and take over 
the whole system after a sufficiently long time. The time evolution of the energy is shown in Fig. 3, confirming our 
previous determination of T m . We therefore conclude that the CTLS has a first order thermodynamic phase transition 
at T m = 1.30. 

Even though the energy of a configuration may be considered as an index of the degree of crystallization of the 
system, it is certainly not a very sensitive one. Given that we know the crystal in this system, and thanks to its low 
degeneracy, we introduce some methods to directly quantify the amount of crystalline order in our samples. This can 
be done either by measuring the typical crystal droplet size £, or by measuring the total crystal mass m. To measure 
£ we observe that both crystals have a periodic structure with a unit cell of size 5. Thus, the Fourier transform of the 
spatial correlation function, 

= do, 

has a peak at k = 2ir/5. The width A of this peak is proportional to the inverse of the typical crystalline domain 
size. Wc want to extract from A a measure of how large is the typical crystal domain, a question which is better 
answered by comparing it to the system size. So we actually define £ as a density, i.e. £ = 1/(AL) G [0 : 1]. It is 
important to note that this definition of £ has a limited power of resolution: in a random configuration wc find on 
average that £o = 0.05, which can therefore be considered as the effective zero for £. When we think of a configuration 
entirely made up of tiny mismatched crystallites we realize that £ is not enough as a measure of crystal order: in this 
case £ would be very small, even though the total amount of crystal would be large. The crystalline mass to can be 
measured by counting the fraction of crystallized spins. To this aim, we have to define an elementary unit, which is 
large enough to have a small probability to be formed randomly, but small enough to be sensitive to small amounts 
of crystal. We choose the nine-spin unit shown in Fig. 1 (right): in the pure crystal each minus spin is surrounded by 
eight pluses, and in the unit cell we have 5 of these elements (which clearly overlap). To count all the spins in the 
cell we must multiply the number of these units by 5. Thus, we define the crystal mass density m G [0 : 1] as the 
number of these elementary units times 5 normalized by L 2 . As for £, we have a random average value of the mass, 
mo = 0.01, which is the effective zero of to. Had we chosen the 5x5 unit cell, the accuracy would be higher, but the 
sensitivity would decrease as well. In order to test the definition of £ and to, we measure their time evolution in the 
experiment described above for the determination of T m . This results are shown in Fig. 4 and Fig. 5. Note that in 
the liquid phase at T = T m + 0.01 we have to ~ 0.15 ^> m (whereas £ is compatible with zero), meaning that also 
in the liquid above T m there is a certain amount of short range order. From this experiment we see that both these 
quantities are excellent configurational indicators of the amount of crystalline order in the system. 

We finish this section by noting that if we extrapolate Elq(T) and Flq(T) to low temperatures, we can find the 
temperature T s where the entropy of the supercooled liquid becomes equal to the entropy of the crystal (which, in 
this temperature range, is to a very good degree of approximation equal to zero). We find 

T s = 0.91. (11) 

Therefore, T s is the temperature where the Kauzmann paradox occurs [10]: below T s the supercooled liquid entropy 
would become smaller than that of the crystal, that is, quite unacceptably, smaller than zero. This temperature is 
usually called Kauzmann temperature in the literature, a denomination we find somewhat unsatisfactory given that 
Kauzmann himself did not believe the extrapolation was really meaningful. Indeed, as we shall see, the Kauzmann 
paradox is avoided in the CTLS. 
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FIG. 3. Energy per spin as a function of time below and above T m . L = 100. 




time 

FIG. 4. Domain size £ as a function of time below and above T m . Same temperatures and symbols as in Fig. 3. 
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FIG. 5. Crystal mass m as a function of time below and above T m . Same temperatures and symbols as in Fig. 3. 



IV. DYNAMICS 

The aim of this section is to study the dynamical behaviour of the CTLS in the supercooled liquid phase (i.e. for 
T < T m ) and in particular to analyze how the relaxation time increases when decreasing the temperature. As already 
anticipated in the Introduction, there is a marked difference in the behaviour of the system if our observational 



6 



time window is longer or shorter than a given time t sp ~ 10 4 MCS. For t < t sp the supercooled liquid is stable 
at all temperatures where equilibration may be achieved: in this time regime we can either study the equilibrium 
supercooled liquid or the off-equilibrium glass, formed by a sufficiently fast cooling. On the other hand, if we are able 
to observe the dynamics on time scales t > t sp we discover that the supercooled liquid loses stability at low enough 
temperatures, and the system enters a completely different regime. Therefore, it is useful to divide our dynamical 
investigation into two parts: first we shall assume that our maximum experimental time t CKp is shorter than t sp and 
we will measure the equilibrium properties of the supercooled liquid. Second, we will consider times larger than t sp , to 
analyze the loss of stability of the supercooled liquid. The scenario that will emerge from this analysis is the following. 
Our system, if analyzed on short time scales exhibits all the features of typical fragile structural glasses, posing the 
usual questions about the nature of the glassy state at low temperature. Only when looking at longer time-scales it 
is possible to unveil the behaviour of the system and give explanations which would otherwise be out of reach. 



A. Moderate times dynamics: t < t. 



In this sub-section we will assume that the longest experimental time available to us is i exp < t sp ~ 10 4 MCS. To 
probe the dynamics of the CTLS we measure the equilibrium relaxation time r as a function of the temperature T. 
This can be obtained from the decay rate of the (normalized) equilibrium spin-spin correlation function, defined as, 

rut \- *(*™+*))-( 8 ) 2 n9 x 
<^\i,t w ) — j - ! 

where (• • •) indicate averages over spins and initial configurations. We quench the system from infinite temperature 
down to the target temperature T and wait for it to equilibrate. At equilibrium the dependence on t w disappears and 
C = C(t). This is our first equilibration test. In order to extract a relaxation time from the correlation function we 
fit it with a stretched exponential, which is the expected equilibrium behaviour for C(t) in super-cooled liquids 

C(t) = exp [-(t/r)"] . (13) 

In Fig. 6 we plot C(t) for three different temperatures, together with their stretched exponential fit. To be sure 
that the system has equilibrated we check that at each temperature the equilibration time that we have waited is 
much larger than the relaxation time r, as obtained from the fits (13). More precisely, since C{t) decays to zero in 
approximately 20t we require t eq > 20r. 




The relaxation time as a function of the temperature is plotted in Fig. 7. We note that there is no dynamic signature 
of the melting transition T m . Below T m the supercooled liquid is metastable with respect to the crystal, and there 
is a sharp increase of the relaxation time. In the inset of Fig. 7 we can see that the Kohlrausch [3 exponent of the 
stretched exponential fit decreases with the temperature, as it happens for realistic models of liquids (see, for example, 
the data of [12]). 



7 




In the temperature regime of Fig. 7 the relaxation time can be fitted very well by a power law, 

with T c = 1.06, 7 = 2.29. 



(T - T c y. 



(14) 



This fact suggests that the CTLS is a fragile system [13]. To support this conclusion, we compare the CTLS with 
an Arrhenius/strong system. In [14] it has been noted that the PQ model in d = 2 is such a system, and we 
therefore compare the two models by plotting the logarithm of r as a function of T*/T, where T* is an normalization 
temperature, arbitrarily defined as the point where r(T*) = 1000. 



1000 




FIG. 8. Angell plot for the CTLS and the PQ model in 2d. 



This Angell plot [13] is shown in Fig. 8: it is clear that the relaxation time of the 2d-PQ model time is very well 
described by an Arrhenius form, r ~ exp(A/T), typical of strong systems, while the CTLS displays the curved shape 
of fragile glasses. To give a quantitative measure of fragility we need to fit the relaxation time in a Vogel-Fulchner- 
Tamman form, 

The fit works pretty well for the CTLS and gives T = 0.76 for T E [1.2 : 1.7], and T Q = 0.90 for T 6 [1.2 : 1.4]. The 
quantity Kvft — T /A is a measure of the degree of fragility of the system [13]. For the two models we find, 



8 



2dPQ : KypT = 0.017, 
CTLS : K VFT = 0.095, 



(16) 
(17) 



consistent with a high fragility of the CTLS. In [6] it has been shown that in d = 3 the PQ model is very well fitted 
by a power law, suggesting that it is a fragile model as the CTLS. Incidentally, we note that an Arrhcnius fit of the 
2rf-PQ correlation time gives a barrier size A = 5.8, at variance with the results of [14], A = 2. Unfortunately, we 
were not able to explain this difference. 

From Fig. 6 we note that the correlation function C(t) relaxes to zero in a time of order 20r, and this, as we 
have seen, is also the order of magnitude of our equilibration time. Thus, for T < 1.25 we are no longer able to 
equilibrate the system within our experimental time t cxp = 10 4 , therefore this is our glass transition temperature 
within this time regime. If, however we give up the requirement to be at equilibrium, we can of course explore 
the low temperature region. For quenches below T = 1.25 the system should be considered, according to usual 
definitions, as an out of equilibrium glass. In fact, out-of-cquilibrium glasses can be formed when cooling the liquid 
at low temperatures at a rate fast enough to prevent equilibration (the quench from infinite temperature is just the 
particular cooling protocol with infinite cooling rate). In Fig. 9 we show the results of cooling experiments with a 
linear cooling schedule T = T in — rt with cooling rate r = dT/dt and T in — 1.67. We plot the system energy as a 
function of T for different cooling rates. In the same graph we report the equilibrium supercooled energy Elq(T), 
which of course must be extrapolated below the last equilibrium temperature. As expected, at high temperatures 
even for relatively fast coolings the system stays in equilibrium, whereas E(T) departs from the equilibrium curve 
at lower temperatures. The slower the cooling, the lower the temperature where the energy freezes in the glassy 
phase, and this same temperature is the glass transition T g (r) corresponding to that particular cooling rate. As in 
the measurement of r, we have a limit to our experiment set by our maximum experimental time t c ^ p : to bring the 
system through the interesting temperature region of width T G [0.8 : 1.3], within a time t cxp we can afford at most 
cooling rates of the order of 5 x 10~ 5 . 




The energy of the configurations reached at T = upon cooling, as a function of the cooling rate, is a quantity 
which is considered as an important indicator of glassy behaviour. What we find for r e [5 x 10~ 4 , 0.01] is a slow 
dependence of E — Eqs on r, compatible with a power law with low exponent or with logarithmic behaviour (see 
Fig. 21 and section V for a wider discussion of this point), which is very similar to what observed in real glasses. 

Summarizing the results of this subsection, what we observe in the CTLS for times smaller than t sp ~ 10 4 MCS is a 
phenomenology typical of fragile structural glasses: in the supercooled liquid phase, T < T m , the correlation function 
displays stretched exponential relaxation. The relaxation time r grows very sharply and is well fitted by a power law. 
Moreover, the exponent (3 of the stretched exponential fit decreases with decreasing temperature. A VFT fit of the 
relaxation time gives a fragility index much larger than a typical Arrhenius/strong spin system. Finally, on rapid 
cooling, the system goes out of equilibrium remaining stuck in the glassy phase, at a glass transition temperature 
which is lower the smaller the cooling rate. 
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B. Long times dynamics: t > t sp and the metastability limit T ap 

We analyze now the behaviour of the system for times longer than t sp = 10 4 . In Fig. 10 we plot the energy as a 
function of time for three different temperatures below T m , with random initial condition. 

At the higher temperature, T = 1.26, the system relaxes in the supercooled liquid and remains in this phase up 
to our actual experimental time, t cxp = 2 x 10 6 . On the other hand, for T = 1.23, after initial equilibration in the 
supercooled liquid phase, the system makes the transition to the crystal. This fact can be appreciated in Fig. 11, where 
we show the behaviour of the crystal domain size £, which is clearly highly correlated with the energy. Note that at 
T = 1.23 the supercooled liquid survives up to 10 5 MCS when the transition to the crystal starts, and crystallization 
is completed in 10 6 MCS. At T = 1.18, on the other hand, we see that the supercooled liquid lasts roughly 10 4 
MCS, while the approach to the crystal takes more than 10 6 MCS. In other words, the lower the temperature, the 
faster is the departure from the liquid, but the slower is the completion of the crystal. If we further decrease the 
temperature, the plateau corresponding to the supercooled liquid disappears, and the energy decreases steadily toward 
the crystal ground state. The dynamics becomes ever slower in this temperature regime, and the system remains out 
of equilibrium up to very long times (see sect. V). 
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FIG. 10. Energy per spin close to the spinodal as a function of time. The dotted line marks the energy of the crystal. 
L — 100. Error bars are only shown if larger than symbols size. 
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FIG. 11. Domain size close to the spinodal as a function of time. Data and symbols as in Fig. 10. 

These results can be interpreted in terms of nucleation theory [15]. Even in the super-cooled liquid phase the 
system nucleates droplets of crystal with finite probability. Most of these droplets are too small to be stable and 
soon disappear, however a critical size exists such that whenever a droplet of that size is created it will not shrink 
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back, but rather expand and bring the system into the stable crystal phase. This phenomenon, critical nucleation, is 
characteristic of metastability and determines the average life-time of the metastable phase. The time needed by the 
system to nucleate critical droplets, i.e. the nucleation time r nuc , depends on temperature and may be very long. It 
is often not accessible experimentally, so that the loss of stability of the metastable phase cannot be observed. This 
is not the case of the CTLS model: for T = 1.23, r nuc is indeed low enough to detect the loss of stability of the 
super-cooled liquid. Once critical nucleation takes place, critical crystal droplets start to grow until all the system has 
a crystalline structure. However, the growth of the crystal may become very slow if kinetic constraints are present. 
What we observe in the CTLS is that the lower the temperature, the faster is crystal nucleation, but the slower is 
crystal growth. 

We can formalize the loss of stability of the liquid by estimating the nucleation time r nuc . In two dimensions, the 
free energy difference due to formation of a crystal droplet of size £ in the liquid phase is given by [15], 

A(0 = Aa£-BS.fe, (18) 

where A and B are geometric factors depending on the shape of the droplet, <r(T) is the surface tension, and 
8f(T) = Flq(T) — Fcr(T) is the bulk free energy difference between supercooled liquid and crystal below T m . 
The function A(£) has a maximum at a value of the droplet size £*, with an associated value A* of the excess free 
energy: 

4 1 ' 2B Sf(T) ' 1 ' AB Sf(T) ' 1 j 

For £ < the system still needs energy to enlarge the droplet which therefore naturally shrinks its size to zero; for 
£ > on the other hand it becomes favorable for the droplet to grow. Therefore, £* is the critical droplet size and 
A* the free energy barrier to crystal nucleation. The Arrhenius estimate for the time T nuc of crystal nucleation is 
then 

/A 2 a(T) 2 \ 

WT)=r exp^-^L_j. (20) 

A necessary condition for the existence of the supercooled liquid phase at temperature T is that the nucleation 
time Tnuc(T') is much longer than the time to make an equilibrium measurement, i.e. the equilibration time r eq (T). 
In other words, only if T nuc (T) > T oq (T) can we speak of a metastable supercooled liquid. We want now to test this 
relation in the CTLS. To find r nuc we need the surface tension a, which we estimate at a reference temperature T re f 
and assume that depends weakly on T close to T re f. To do this we invert formula (20), 

a 2 = 4BT ref 5 f {T ref ) log[r nuc (T re/ )/r ] /A 2 , (21) 
and after setting a — T re f 5f(T re f)/TSf(T), we can write 

r nuc (T) « r expja log [ T -(?W) j j (22) 

— T nuc(T r ef) a T Q (23) 
» Tnuc(T ref ) a (24) 

where the last approximation is valid when T is near T re f , since then awl. To have a result as accurate as possible 
close to the loss of stability of the liquid, we take T re f — 1.234, where T nuc = 5 x 10 4 . 

Regarding the estimate of r oq (T), we have already noticed that the correlation function C(t) drops to zero after 
roughly 20 relaxation times, and this we take as a reasonable estimate of the equilibration time (it is likely that we 
are underestimating r cq in this way). The relaxation time is given by the power law fit, eq. 14. 

In Fig. 12 we plot T nuc and r cq vs. T from the above estimates. We see that the nucleation time drops below the 
equilibration time at a temperature T sp — 1.22, which therefore marks the metastability limit of the supercooled 
liquid phase. Below T sp the liquid can only be observed on time scales shorter than the equilibration time, meaning 
that what we are actually observing is an out-of-cquilibrium glass. Therefore, the equation 

^nuc (?sp) = T C q(T S p) (25) 

can be considered as a definition of the effective spinodal temperature T sp for the metastable liquid phase. Note that 
the two curves cross at a time t sp ~ 10 4 MCS: this is the reason why if our maximum experimental time is smaller 
than t sp we do not detect this metastability limit. As we have already said, the glass transition temperature T g is 
fixed by the experimental time t cxp via the relation r cq (T g ) = i cxp . Therefore, if T sp < T g the metastability limit 
cannot be observed, whereas we detect it if T g < T sp . 
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FIG. 12. Nucleation and equilibration times vs temperature. 



The presence of a lower metastability limit T sp is clearly visible also in cooling experiments. If our experimental 



time exceeds t sp = 10 4 MCS, we can afford smaller cooling rates than 



T e [0.8 : 1.3]. The results are reported in Fig. 13. For r > 



sp 



AT /t sp to explore the relevant region 



sp 



5 x 10 the energy departs from the equilibrium 



line remaining above it, while for r < r sp it drops below it. We can understand qualitatively this behaviour in terms 
of equilibration and nucleation processes. For r > r sp , given the fast cooling rate, the system has not the time to 
equilibrate at low temperatures where the equilibration time becomes very large and its energy therefore remains 
higher than the equilibrium one. On the contrary, for r < r sp the cooling is so slow that at low enough temperatures 
the system has the time to nucleate critical crystal droplets and its energy is driven toward the ground state one, thus 
becoming smaller than the supercooled liquid energy. In both cases, however, the system eventually remains stuck in 
an off equilibrium phase. 
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FIG. 13. Cooling experiments above and below the spinodal cooling rate r sp . Cooling procedure as in Fig. 9. L = 100. 



This experiment poses an interesting question about the actual nature of the out-of-cquilibrium glassy phase. For 
fast cooling rates the system remains stuck at a high energy level, into a clearly disordered, or glassy-like, phase 
(Fig. 14, upper left). On the other hand, for the slowest cooling rates the system approaches the ordered crystal 
phase, even though it remains stuck in an out-of-equilibrium polycrystalline configuration, as Fig. 14 (lower right) 
shows. The decrease in energy, and the associated increase in crystalline order that we observe as we slow down the 
cooling is a continuous process, making it very difficult to distinguish what is a glass and what is a highly disordered 
polycrystal: every possible structural way to detect crystalline order in the upper left configurations of Fig. 14 would 
give the same result as for the liquid phase. Note, finally, that all the configurations in Fig. 14 are stable local minima 
of the Hamiltonian reached at T = 0. 
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FIG. 14. Local minima of the CTLS reached at T = with different cooling rates, from r = 2 x 10 3 (upper left), to 
r = 2 x 1CT 6 (lower right). L = 100. 



This structural continuity from glass to polycrystal can be appreciated more quantitatively in Fig. 15, where we 
plot the energy of the asymptotic configurations the system reaches at T = 0, as a function of the corresponding 
crystalline domain size £ and crystal mass m, for various cooling rates r. All the points on this plot corresponds to 
configurations obtained with a well-defined cooling experiment. Moreover, all these configurations are (local) minima 
of the energy. As we can see, the data interpolate without gaps between crystalline (or polycrystalline) configurations 
with low E and high £ or m, to glassy configurations with high E and low £ or m. In this sense we may say that in 
the present system there is no qualitative difference between the off-equilibrium disordered, i.e. glassy, regime, and 
the crystal growth regime. On the other hand, thanks to the presence of the pseudo-spinodal transition, there is a 
clear distinction between equilibrium liquid and crystal phases. 
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C. Metastability limit and relaxation time divergence: T sp vs T c 

It has been shown in [6] that in the PQ model the temperature T c where a power law fit of the relaxation time locates 
a divergence, is very close to the metastability limit T sp . The interpretation of this fact given in [6] is very clear: in 
mean-field the PQ model has a sharp spinodal temperature at the point where the liquid metastable state disappears, 
and this transition is accompanied by a divergence of the correlation times. Therefore, in [6] it is argued that in the 
finite-dimensional version of the PQ model the temperature T c is actually a remnant of the mean-field spinodal. The 
same view has been advanced for off-lattice models: in [16] the behaviour of a supercooled two-component Lennard 
Jones system in d = 2 has been interpreted as due to the influence of a pseudospinodal, while in [17] arguments are 
given to explain why in d = 3 spinodal effects may be much harder to detect. 

The situation is very similar in the CTLS, even though our definition of T sp makes it indispensable to distinguish 
between the two temperatures. The temperature T c where the relaxation time would diverge is not actually a true 
critical temperature: when the temperature is decreased the equilibration time increases and finally becomes as large 
as the nucleation time at T sp , when r cq = r nuc is still finite. Therefore we have by definition T sp > T c . This is simply 
due to the finite-dimensional nature of the CTLS. In mean- field systems a metastable state has infinite life-time and 
is separated from the stable phase by an infinite free energy barrier: the spinodal temperature, which marks the 
disappearance of the metastable state, corresponds to the temperature where the barrier goes to zero, and is naturally 
associated with a divergent time-scale. In a finite dimensional system like the CTLS, on the contrary, the metastable 
phase (the liquid) is separated from the stable one (the crystal) by a finite free energy barrier — see eq. (19) — and 
always has a finite lifetime. An effective spinodal temperature can only be defined by comparing the different and 
finite time-scales present in the system, as done in the previous subsection. However, the physical interpretation of 
[6] may remain valid, i.e. the increase of the relaxation time close to T sp may be an effect of approaching incoming 
instability of the supercooled liquid. 

We stress again that had our maximum experimental time been shorter than t sp = 10 4 MCS, we would not have 
detected the metastability limit at T sp , while we still would fit the equilibrium relaxation time as a power law with 
critical temperature T c . In that case, it would have been harder to realize the connection between divergence of the 
relaxation time and loss of stability. 

D. Kauzmann's resolution of the Kauzmann paradox 

In 1948, W. Kauzmann outlined in a famous paper [10] the paradox now named after him. In the supercooled phase 
the viscosity of the system increases with decreasing temperature, until at the glass transition T g relaxation time 
becomes too long and equilibration cannot be achieved within experimentally accessible times. However, Kauzmann 
noted that if the entropy of a supercooled liquid is extrapolated below T g , it becomes equal to the crystal entropy 
at a temperature T s > 0, and even smaller than zero if extrapolated further. This entropy crisis is never actually 
observed, because the glass transition intervenes before. However, Kauzmann found it paradoxical that it was just a 
kinetic phenomenon (the glass transition) that saved the liquid from a thermodynamic nonsense. 

In the context of the Adams-Gibbs-DiMarzio theory [18], the entropy crisis has however an interpretation: the 
entropy difference between crystal and liquid is related to the configurational entropy S, that is the entropic contribu- 
tion due to the presence of an exponentially high number of different glassy minima. The vanishing of £ at T s signals 
a thermodynamic transition to a new phase, characterized by a sub-exponential number of glassy states, separated 
by infinite free-energy barriers. This picture is exact for some mean-field spin-glass systems [19], and it may be the 
correct resolution of the Kauzmann paradox even for real structural glasses. According to Adam-Gibbs-DiMarzio, 
the entropy crisis is accompanied by a divergence of the relaxation time at T s given by r ~ exp(^4/TX). Expanding 
linearly the configurational entropy close to T s a Vogel-Fulchner-Tamman behaviour is obtained (see eq. 15) with 
T =T S . T and T s have been observed to be quite close in various systems [1], a fact that has been advocated as an 
indirect argument in favour of this interpretation. 

Despite recent analytic and numerical work supporting the entropy crisis scenario [21,20], there is another way to 
avoid the Kauzmann paradox, which, interestingly enough, was proposed by Kauzmann himself [10]. He rejected the 
idea of a thermodynamic glassy phase, and of a transition at T s . What Kauzmann hypothesized is the existence 
of a metastability limit of the supercooled liquid phase, below which crystal nucleation becomes faster than liquid 
equilibration. More specifically, he defined an effective spinodal temperature T sp > T s below which "the free energy 
barrier to crystal nucleation becomes reduced to the same height as the barrier to simpler motions". Below T sp the 
supercooled liquid is operationally meaningless and thus the paradox is avoided. 

What we have shown in this paper is that the CTLS is a system with a super-cooled phase and a fragile-glass 
phenomenology that precisely reproduces the scenario hypothesized by Kauzmann. If we look at the behaviour of the 
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extrapolated excess entropy Slq — Sen we find that in this model T s and To are actually very similar (see sections 
III and IV- A). However, despite this fact, no entropy crisis does actually take place. The metastable liquid ceases to 
exist below T sp = 1.22 and thus extrapolations of thermodynamic quantities in the supercooled phase below T sp are 
meaningless. Therefore, the paradox is avoided: nucleation takes over when the relaxation time becomes larger than 
the nucleation time. Of course, we can still speak of an out of equilibrium glass below T sp , which can be obtained if 
we cool fast enough our system. 

As we have already stressed, the metastability limit may prove impossible to observe experimentally if the equili- 
bration time at T sp is much larger than the experimental time, that is if T sp <C T g . One is then tempted to suggest 
that the Kauzmann's resolution of the Kauzmann paradox that we have explicitly analyzed in the CTLS model pos- 
sibly applies to many other systems, but is actually out of experimental recognition. However, we must be careful in 
judging how general this behaviour can be. The rest of this Section is devoted to a discussion of this point. 

1 e+08 n— — i — ' — i — ■ — i — ' — i — ' — i — ' — i — 



gj le+06 



c 

% 10000 
a 
o 

"o 

3 

c 100 



I I I I I I I I I I I I L_ 

0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 

temperature 

FIG. 16. Extrapolated nucleation time of the CTLS model vs. temperature. A minimum is clearly visible at T — 0.98. 



In Fig. 16 we see that the nucleation time r nuc of the CTLS has a minimum at a temperature T = 0.98, below 
which it starts increasing again. This is simply due to the balance between the free energy barrier A*(T) and T in 
the Arrhenius factor, and it is certainly a very model-dependent property. 1 In the CTLS all this is of little relevance, 
because T nU c and r eq cross at a much higher temperature, as shown in Fig. 12 and also in Fig 17-Left. Yet, in general 
the interplay between a(T), Sf(T) and T may be such that a minimum of r nuc (T) actually occurs and that at this 
minimum r nU c S> r cq . However, even in this case it may still happen that the curves cross below the minimum of T nuc , 
as shown in Fig. 17-Centre. In this case the Kauzmann paradox would still be resolved by the metastability limit, 
exactly as in CTLS. Note that this particular scenario may actually prove to be quite deceiving: at temperatures 
smaller than the minimum of the nucleation curve one can think that the 'dangerous' zone for crystal nucleation has 
been safely crossed and that crystallization is therefore ruled out at lower temperatures. When T g is met at lower 
temperatures, one would need some very careful extrapolations of both r nuc and r cq to recognize that a metastability 
limit actually exists at T sp <C T g . 

But is it possible that T nuc 3> r cq at all temperatures? The answer is certainly yes for strong liquids, where T s = 0, 
and thus r eq stays finite for T > 0; but this case is of little interest for us since then the Kauzmann paradox does not 
arise. In fragile liquids, however, where r eq is supposed to diverge at T s , the condition for the existence of the liquid, 
Tnuc 3> T eq , requires that r nU c diverges at a temperature T nuc >T S . In our opinion it is possible that T nuc = T s , since 
at low temperatures the contribution of t in eq. 20 can be important due to kinetic constraints, as shown in [22]. 
In nucleation theory, this prefactor is taken as proportional to the viscosity, in which case r nuc would diverge when 
the viscosity does, that is at T s . But it is hard to imagine how T nuc could diverge strictly above the thermodynamic 
transition temperature T s . Therefore, it seems to us that the best one can have is T sp = T s , as in Fig. 17-Right. 

Let us stress once again the fact that crystal nucleation is a different process from crystal growth: saying that 
nucleation is fast, is not the same as saying that the critical crystal droplets will grow quickly, since the speed of 



1 We have assumed a constant surface tension close to T re f — 1.23, and thus we cannot trust too much r nuc far from that 
point. Moreover, we have ignored any possible temperature dependence in the prefactor of Arrhenius's law (20). 
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growth will strongly depend on kinetic considerations. From an experimental point of view, it is clear that having a 
very long nucleation time (longer than the experimental time) and fast growth, is not much different from having a 
fast nucleation of very small, undetectable, crystal droplets, which however grow exceedingly slowly. In order to avoid 
crystallization and form the glass we have either to avoid nucleation, or to ensure that crystal growth is kinetically 
blocked. To clarify this point, we can ask how long is the time t x the system takes to develop a substantial amount 
of crystalline order, large enough to prove that crystallization has started. We can answer this question by following 
the time evolution of the crystal mass m(t) at various temperatures, and defining t x from the relation m(t x ) > mth, 
where m t h is an arbitrary threshold. The function t x (T) is normally called in the context of supercooled liquids Time- 
Temperature- Transformation (TTT) curve, and it is the only experimentally available means to check how likely is 
crystal formation at a given temperature. Note that the value of t x (T) will strongly depend on m t h, and it cannot 
therefore be quantitatively compared to r eq . 
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FIG. 17. Pictorial view of three different scenarios: Left: T nuc and r eq cross at a spinodal temperature T sp , which is larger 
than the temperature where r nU c has a minimum, as in the CTLS. Centre: T nuc and r eq cross at a spinodal temperature T sp , 
which is smaller than the temperature where r nuc has a minimum. Right: r nuc and r eq both diverge at the same temperature 
T s — T sp . The circles indicate the spinodal temperature T sp . 
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FIG. 18. TTT curves for two different values of the threshold value m th . L=100. 

In Fig. 18 we plot the TTT curves of the CTLS for two different values of m t h- As in realistic supercooled liquids, 
the TTT curves have a minimum at a temperature T x : here crystallization has the fastest rate, in the sense that at 
T x we have to wait the minimum possible time to reach a level m t h of crystallization. The origin of this minimum is 
obvious in CTLS: for high temperatures t x is large because crystal nucleation is very slow, i.e. because T nuc is very 
large, whereas for small temperatures t x is large because crystal growth is very slow, which is essentially due to the 
high viscosity of the system (i.e. to the same reason why, if the liquid still existed, r cq would be very large). Thus, 
in CTLS the two branches of the TTT curve correspond to two different phenomena, that is nucleation and growth. 
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Not surprisingly we find T x ~ T sp . On the other hand, in a system where the scenario depicted in Fig. 17-Centre 
holds, also the left branch of the TTT curve for T < T x would be due to slow nucleation [23,24], until T sp < T x is 
reached. In this case the minimum of the TTT curve would be an effect of the minimum of r nuc (T) described above. 
Below T sp , however, the TTT curve would again increase due to slow crystal growth. 

It may be argued that there is little experimental difference between the case where crystallization at low T is blocked 
because of arrested nucleation or arrested growth. However, from a theoretical point of view these two situations are 
quite different: in the first case the equilibrium liquid phase exists, at least on certain time scales, whereas in the 
second case it does not, the only liquid phase is the out-of-equilibrium glass and the Kauzmann paradox is avoided. 
We can characterize these two scenarios as equilibrium arrested nucleation vs off- equilibrium arrested growth. 

V. OFF-EQUILIBRIUM BEHAVIOUR: DOMAIN GROWTH AND AGING 

In this section we study the dynamical behaviour of the CTLS after a quench to T < T sp . We shall first focus 
on one-time quantities, and in particular on the behaviour of energy, crystal domain size and crystal mass. Then 
we will analyze two-time quantities to understand whether our system exhibits aging, and of what kind. The key 
result of this section is that there are two time regions for off-equilibrium dynamics. In the earlier time regime there 
is the formation of a large number of critical crystal droplets all over the system. As we have seen, this mixture of 
crystallites cannot be structurally distinguished from a glass or liquid configuration, as long as the crystallites are 
very small. Their growth is controlled by activation, and thus this regime will last longer the lower the temperature. 
In this phase, which we call bubbling, off-equilibrium dynamics is similar to what found in structural glasses. The later 
time regime kicks in when the energy barriers to crystal growth have been crossed, so that crystal domains are large 
enough to be structurally recognizable. This phase we call coarsening, and off-equilibrium dynamics in this regime is 
consistent with models where simple coarsening occurs, as the Ising model. Not surprisingly, the minimum time we 
have to wait to see the coarsening regime in our system (and thus realize that crystal domains are actually growing) is 
as long as the time needed to directly see the loss of stability of the liquid, that is t sp . In other words, the dynamical 
checks for the growth of the crystal are as inefficient as the structural ones we considered in the former Section. 

A. One time quantities 

As seen in the previous sections, when quenched to temperatures in the range T sp < T < T m , the system relaxes 
into the supercooled liquid phase and here remains up to times of order r nuc (see Fig. 3). On the other hand, below 
the spinodal temperature T sp the supercooled liquid no longer exists and after a quench the energy very slowly evolves 
toward its crystalline ground state value. The rate of decay of the energy gives a useful characterization of the off- 
equilibrium dynamics. In Fig. 19 we report the behaviour of the excess energy SE = E — E GS as a function of time, 
together with power law fits of the form 5E = At~ v . For all temperatures T < T sp we find an initial time regime 
where v < 0.2, which is compatible with what found in some model systems of structural glasses (see for ex. [25,26]). 
This first time regime is followed by a second late time regime where the energy decreases faster, with v ~ 0.5. 
However, the data show that the crossover to this second time regime takes place later the lower the temperature. At 
T = 0.67 we are unable to enter this late time region within our maximum experimental time, and we only observe 
the slow regime. 

The existence of these two different regimes can be related to the presence of activated coarsening, as is commonly 
found in frustrated lattice spin models. [4,5]. Below the spinodal temperature the system rapidly nucleates small stable 
crystal domains, which tend to grow. However, due to the high viscosity, local spin rearrangements arc difficult, and 
moving a domain boundary may be rather costly. Thus, contrary to the case of standard coarsening in the Ising 
model, when frustration is present domain boundaries are pinned, and domain growth requires overcoming some 
energy barriers. The system therefore needs activation to grow crystal domains, and this may lead to a non-standard 
decay of the energy as a function of time. 

More precisely, according to [27] we must distinguish two classes of activated coarsening, depending on whether 
or not the maximum barrier A to domain growth depends on the linear domain size L£(t), and therefore on t itself. 
In the first class of systems A does not depend on In this case the system evolves slowly within a time scale 

*A ~ exp(/3A), while for logt > log^A barriers are no longer a problem, and a standard power law behaviour is 
recovered, i.e. E — Eqs ~ t -1 / 2 (we are only considering systems with nonconserved order parameter). On the 
contrary, if the barrier to domain growth is proportional to £(t), then it is larger the longer the time. In this second 
class of activated coarsening one expects a logarithmic asymptotic dependence of the energy on time. 
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FIG. 19. Excess energy as a function of time after a quench at T < T Bp . Full lines: power law fit of the early (slow) time 

1.00. Dashed lines: power law fit of the late (fast) time 
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FIG. 21. Energy of the configurations reached at T = during cooling experiments (see Fig. 13) as a function of the cooling 
rate, together with power law and logarithmic fits. L = 100. 
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The data regarding the energy in Fig. 19 show that, after an early slow regime, the standard i -1 / 2 decay of the 
energy is recovered. This fact suggests that our model belongs to the first class of activated coarsening, and that 
the early time regime is indeed connected by activation to a constant energy barrier to domain growth 2 . A similar 
conclusion can be inferred from the behaviour of In figure 20 we show for three temperatures below T sp the 

behaviour of £ as a function of time. Also in this case for the higher temperatures we find quite clearly two different 
regimes: a first one in which £ does not increase, but rather fluctuates around its effective zero, and a second one 
where it increases quite sharply. The statistics is not sufficient to envisage a clear asymptotic dependence, however a 
power law fit gives for the late time regime an exponent close to 1/2. 

The behaviour of the system when we perform a cooling experiment is consistent with the interpretation above. 
The configurations reached at T = for low cooling rates exhibit the typical pattern of a coarsening system (Fig. 14) 
and are stable, indicating that coarsening is indeed activated. Besides, if we plot the excess energy SE as a function 
of the cooling rate r, once again we find two regimes (fig. 21): for fast coolings (short time) we have a slow energy 
decay, which can be fitted either with a logarithm, or with a power law 5E ~ r", with v ~ 0.15. On the other 
hand, for slow coolings (long time) we find v ~ 1/2, as in standard coarsening for the Ising model. This behaviour of 
SE(r) is what we expect for activated coarsening of the first class and can be qualitatively explained in the following 
way. As we have seen, below T sp the system exhibits a non trivial dynamical behaviour where the energy evolution 
is determined by the formation and growth of pinned crystal domains. In a cooling experiment it is then crucial how 
much time the system spends below the spinodal temperature before remaining trapped in the asymptotic state. If we 
call Tf the temperature where the system reaches the asymptotic T = energy value 3 , then at a given cooling rate 
r the time elapsed from T sp to Tf is given by t(r) = (Tf — T sp )/r. During this time the energy of the system evolves 
approximately as described previously in this section, i.e. 5E = At~ v . Therefore we have E(r) — E GS ~ t(r) _,y ~ r v 
and the two regimes observed in Fig. 21 are nothing else than a manifestation of the two regimes already discussed 
for Fig. 19. The exponent v corresponds to the first regime value if t < t sp , that is if r > r sp , while for r < r sp the 
standard value v = 1/2 is recovered. 

B. The bubbling and coarsening regimes 

The two regimes we have presented above can be explicitly related to different patterns of domain dynamics. In 
order to do this we compare the different time dependence of the domain size £ and crystal mass m 4 . Up to now we 
used both quantities as equivalent markers of crystal order, but in fact in the growth phase their role is quite different. 

After a quench to T -C T sp we expect that fast nucleation will lead to an increase in m due to the formation of 
many crystallites, but that the typical domain size £ will remain very small, actually smaller than our resolution. This 
picture is confirmed by Fig. 22, where we plot £ and to as a function of time for two different temperatures below 
T sp . In both cases we see that the initial increase of m is faster than that of £. For T = 1.0, which is quite close 
to T sp = 1.22, we can see a rise of £ only for times longer than 10 5 MCS, while for T = 0.67 there is practically no 
increase of £ up to t = 10 6 MCS, while a clearly detectable increase of m can be seen. 

To further support this picture we plot in Fig. 23 the mass to as a function of domain size £ for many different 
temperatures below T sp , and for different times. All these data collapse on a master curve, which shows that for short 
times (small £ and to) the crystal mass grows much faster than the domain size, whereas for long times m tends to 
saturate, and £ start increasing rapidly. 

We can therefore define two time regimes for crystal domain growth. For shorter times, in what we may call a 
bubbling regime, many small crystallites form at the expenses of a 'liquid' background, their linear size growing very 
little with time, whereas the total crystal mass grows quite sharply [28] . In this phase domains are so small that they 
are not even well defined, in the sense that their linear size L£ is of the same order as the width a of the interfaces 
(or domain walls) among them. In this regime the excess energy is dominated by the regions not yet crystallized, and 
the interfacial and domain contributions are of the same order. For longer times, the largest part of the system is 
crystallized, so the mass to saturates, and the small domains start growing one at the expenses of each other. 



2 We note that in the models studied in [4,5], a detailed analysis of the barrier to growth showed that these two models actually 
belong to the second class of activated coarsening. 

3 Tf can be estimated as the temperature where thermal activation becomes extremely rare even for simple rearrangements, 
i.e. S/Tf S> 1 where S is the smallest barrier to kinetic motion. 

4 We recall that both these quantities are define as densities. The actual linear length of a domain is L£, and the total Crystal 
mass is L 2 m. 
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FIG. 22. Mass and domain size as a function of time for two different temperatures. L — 100 
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FIG. 23. Crystal mass as a function of the domain size parametrically in t, after a quench at various temperatures. Inset: 
excess energy as a function of the domain size. Full line: l/£ fit. All runs are 2 x 10 6 MCS long, L — 100. 



Therefore a proper coarsening regime starts with a clearly increasing £, and L£ 3> a. Now domains are well defined 
and the excess energy is concentrated in the interfaces. These are now much thinner than the domain size and can 
be considered as lower dimensional manifolds. This must give, 

as in standard coarsening [29,27]. In Fig. 23 we plot the excess energy E(t) — Eqs vs. £ for a number of temperatures 
and times. Consistently with our conclusions, for shorter times there is a drop in energy with basically no increase 
in £. On the other hand, for £ > 0.2 the l/£ law fits reasonably well the data, suggesting that this is the late-time 
coarsening regime. Note that £ > 0.2 implies m > 0.5 (see Fig. 23). We could thus schematically say that the 
coarsening regime sets in when the total crystal mass exceeds 50%. This definition of the two regimes is consistent 
with the behaviour displayed by E(t): for T = 1 the system enters the coarsening regime (i.e. m > 0.5 and £ > 0.2) 
for t ~ 40000 (see Fig. 22), precisely on the same time scale the energy (Fig. 19) passes from a power law behaviour 
with low exponent to a faster t -1 / 2 behaviour. For T = 0.67 on the other hand £ < 0.2 for our longest time, and 
correspondingly the energy always obeys a slow power law decay. 
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C. Two time quantities 



We now focus on the behaviour of two-time quantities, for quenches below T m . In equilibrium dynamics, correlation 
functions C(t,t w ) = (A(t w )A(t)) and their associated response functions R(t,t w ) = 6(A(t)) /5h(t w ), (A is a generic 
variable and h is the field conjugate to it) depend only on the difference t — t w . This is not the case when the system 
is out of equilibrium. The aging regime is that in which one-time quantities (such as energy or density) are stationary 
or change very little, but two-time quantities depend explicitly on both times t w ,t [30]. 

An important quantity which allows to distinguish among different kinds of aging behaviour [30] is the fluctuation- 
dissipation (FD) ratio, 

_ R(t,t w ) 

dc(t,t w )/dt w • y ' 

In equilibrium, the fluctuation-dissipation theorem (FDT) states that X = 1, while out of equilibrium the FDT is 
violated and X depends on t, t w in a non trivial way. It has been conjectured [31] that during aging X depends on 
time only through the correlation function, conjecture that has proved valid in all systems studied to date. If one 
considers the susceptibility, or integrated response, 

X(t,t w )=j R{t,s)ds, (28) 



which is much easier to compute numerically, then assuming X — X[C(t,t w )] one finds from (27) 

r c(t.t) 

'C(t,t„ 



fC(t.t) 

X(t,t w )=(3 X(C)dC. (29) 

Jc(t.t,„) 



Therefore one can plot \ vs - C a ^ fixed t w , and extract X from the slope of the curve. In equilibrium X = 1 and 
the curve is a straight line \(t — t w ) = /3[C(0) — C(t — t w )], with slope — l/fc^T. For a system out of equilibrium one 
finds in general two regimes: a short time quasi-equilibrium regime (t — t w <^t w ), corresponding to equilibration of 
fast degrees of freedom, where the curve follows the equilibrium straight line and FDT holds; and the aging regime at 
larger times (t — t w <^t w ) where x{t,t w ) departs from the equilibrium line. In particular, for structural glass models 
studied so far via numerical simulations [25,32] it has been found that after leaving the equilibrium line, the curves 
follows another straight line, with smaller (in absolute value) slope. 

To study the FD relations in CTLS we follow (as has been done before in other spin systems [33]) the evolution of 
the staggered magnetization (r)i = ±1 are quenched random variables), 

m a (t) = -i J") (30) 

i 

when addi ng a t erm —hNm s to the Hamiltonian. If the field h is applied at time t w , the 77-averaged susceptibility 
X(i; tw) = m s (t)/h is related to the spin-spin correlation function: 

C s (t,t w ) = ^'£2(s i (t w )s i (t)). (31) 

i 

We have used fields h = 0.1 ksT and h = 0.15 fc^T, and checked that they are within the linear regime. The results 
presented here correspond to the highest field. 

We look for aging behavior in the regions T sp < T < T m and T < T sp . As we have seen, in both these regions 
the system has trouble equilibrating, though for different reasons. We first quench the system to T = 1.25, which is 
above the spinodal, and plot \ vs. C s for different waiting times (Fig. 24). Except perhaps for the shortest t w , the 
curve obeys the FDT: we obtain a straight line with slope —l/ksT, independent of t w . No aging is observed, and the 
behavior is compatible with a system in equilibrium, except that we know from thermodynamics that the equilibrium 
phase at this temperature is the crystal: we are witnessing the metastable liquid, which at this temperature lasts much 
longer than its relaxation time, and can thus be observed. 

On the other hand, a quench to T = 0.67 <SC T sp shows aging very clearly, as the correlation function plotted for 
different waiting times shows (Fig. 25). The FD ratio (Fig. 26) also shows very clearly that the system is out of 
equilibrium. In fact, we know that the liquid is unstable at this temperature, and that the crystal phase is growing. 
In a system undergoing coarsening, like the Ising model quenched below T c , the susceptibility rapidly saturates at 
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the value corresponding to the low-temperature phase, because the contribution of the interfaces is negligible. The 
correlation function, however, continues to drop because of the movement of the domain walls. Thus the FD curve is 
flat for small values of the correlation function [34,33,35]. We would expect this same behavior whenever the dynamics 
is given by growing domains with bulk equilibrium properties, separated by sharp interfaces. However, the FD curve 
does not become flat in the CTLS at T = 0.67, even for our longer waiting times. Instead, the system shows FDT 
violations similar to those commonly exhibited by structural glass models [32,25] and certain mean field spin-glasses 
(solved by a one step replica broken solution, 1-step RSB) [31]. There are two regimes: a quasi equilibrium one with 
slope —(3 and an off-equilibrium one with slope — (3 c g = —Xj3. Note that the FD curves hardly depend on t w , at least 
for the times the simulations can reach. 
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FIG. 24. Parametric plot of x(Mu>) vs C(t,t w ), at various values of the waiting time t w for T = 1.25, above the spinodal 
temperature. L — 500. The full line represents the equilibrium relation x = P0- — C). 
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FIG. 25. Correlation function as a function of t, at different values of the waiting time t w . T = 0.67, L — 500. 

Taken at face value, this result would lead to interpret this as an out-of-equilibrium glassy phase (i.e. a liquid 
with weakly- broken ergodicity). To understand it correctly however we need to recall the discussion of the previous 
section. We have seen that at low temperatures (T < T sp ) the crystal growth dynamics proceeds through two 
qualitatively different time regimes. In the first bubbling regime disordered regions separates small crystal domains 
and the intcrfacial contribution to the susceptibility is still very large. The coarsening regime, where the dynamics is 
dominated by the nucleated regions and the interfacial contribution can be neglected, will only be reached later, when 
the crystal regions have grown enough to overwhelm the interfaces. Due to the extremely slow activated dynamics the 
bubbling regime can last for long, and indeed at very low temperatures it is the only one accessible in our simulation 
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times. If this is the case, we do not expect to see in the response-correlation relation the FDT violations typical of 
coarsening systems, since the coarsening regime has not yet been attained. This is exactly what happens for T = 0.67: 
we can check that the crystal mass for the longest time explored in Fig. 26 is m ~ 0.3 (fig. 22- Right), well below 
the value m <~ 0.5 that characterizes the crossover to the coarsening regime. The 1-RSB-like FDT violation observed 
at T = 0.67 thus occurs in the bubbling regime, while if we could wait longer times we would observe a crossover to 
typical a coarsening behaviour. 




FIG. 26. Parametric plot of \{t, t w ) vs C(t, t w ), at various values of the waiting time t w for T = 0.67, well below the spinodal 
temperature. L — 500. The full line represents the equilibrium relation x = /3(1 ~ C) 

Since the duration of the bubbling regime is shorter for higher temperatures, we expect to explicitly observe an 
FDT crossover to the coarsening regime at higher temperatures. From the former Sections we know that at T = 1.00 
the crossover from the bubbling to the coarsening regime takes place at about 10 4 - 10 5 MCS. In Fig. 27 we plot the 
FD curves for T = 1.00, and the data confirm our expectations. For short waiting times, the curves initially follow 
what seems to be the path of a 1-RSB-like violation. This is the bubbling phase, pretending to be a glass. The curves 
depend strongly on waiting time, and for higher t w evolve toward the normal domain-growth curve (two straight lines, 
with slope —0 at short times and at long times), which is nearly reached for t w = 10 6 MCS. 




FIG. 27. Parametric plot of x(Mu>) vs C(t,t w ), at various values of the waiting time t w for T = 1.00, below the spinodal 
temperature. L — 500. The full line represents the equilibrium relation x = P0- ~ C). 
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VI. CONCLUSIONS 



We have presented a lattice model, the deterministic CTLS, which exhibits a phenomenology typical of glass forming 
systems: it has a melting transition and a low temperature crystalline ground state, but it exists as a supercooled 
liquid for a wide range of temperatures; in the metastable supercooled phase the relaxation time increases in a way 
resemblant of typical fragile glasses and upon fast cooling crystallization can be easily avoided to bring the system in 
disordered glassy-like states. 

The feature that makes this model particularly interesting is that it displays a metastability limit which can be 
explicitly observed within the time scales available for numerical simulations. What we have shown is that the typical 
phenomenology of fragile glass forming systems is entirely compatible with the existence of a metastability limit 
of the supercooled phase, even if it is plausible that, contrarily to the CTLS, for many models this limit cannot 
be observed. A crucial point in this regard concerns the length of experimental times under which measurements 
are performed. By definition the supercooled liquid looses stability at T sp , when crystal nucleation becomes faster 
than liquid equilibration. If the equilibration time at T sp is much larger than the experimental time, however, the 
metastability limit cannot be observed. Since the glass temperature T g is the lowest temperature at which the system 
is equilibrated within the experimental time, this situation takes place when T sp < T g . 

If the system has a metastability limit, below the pseudospinodal T sp the dynamical behaviour basically consists 
of very slow crystal growth. This suggests some criteria to detect whether a system has such a metastability limit or 
not, even when T sp is experimentally inaccessible. For example, one would a priori think that there is a qualitative 
difference between a disordered glassy configuration, obtained by quenching a liquid, and a crystalline configuration, 
however rich in defects this is, and however slow crystal growth may be. Also, since off-equilibrium dynamics of 
simple domain growth has very peculiar characteristics one would hope to detect crystal growth by analyzing the 
dynamical behaviour of the system, and in particular the pattern of FDT violations. What we have shown in the 
CTLS is that actually neither of these criteria is sharp enough, and if our experimental time were not long enough to 
explicitly observe the loss of stability of the liquid at T sp , it would be impossible to make statements about it, cither 
by structural or dynamical means. On the other hand, on shorter time scales the system exhibits an off-equilibrium 
behaviour which is typical of structural fragile glasses. Indeed, we have seen that there is a long time regime where 
FDT violations are 1-RSB like, as in most of fragile glass models studied so far, while the proper coarsening regime is 
attained only at much larger times. Also, there is structural continuity between the low-temperature configurations 
reached via fast coolings and those approached with slower cooling rates, which means that the degree of order in the 
system increases gradually: what appears as a disordered glassy like configuration is just the beginning of a continuous 
process that eventually, on much longer time scales, leads to the crystal. The reason for this behaviour is that below 
T sp crystal nucleation is fast, but crystal growth becomes very slow, with many crystal droplets trying to expand 
in a liquid background. In such a situation distinguishing between a truly disordered glass and a mixture of tiny 
mismatched crystallites becomes very hard, and FDT violation is nontrivial. 

An important consequence of the metastability limit scenario in the CTLS is that the Kauzmann paradox is avoided: 
below T sp the liquid does not exist in any reasonable sense, thus extrapolation of the excess entropy is meaningless. 
This was the resolution of the paradox proposed by Kauzmann himself in his paper. While we cannot claim that this 
is the case for real glasses, (and we have discussed qualitative alternative scenarios in terms of crystal nucleation and 
growth), this system shows that Kauzmann's way out of the Kauzmann paradox ought to be considered seriously. 
In the study of glassy systems the crystal phase is often neglected, assuming that if crystallization is avoided at T m , 
then the crystal does not play any role in the low temperature physics. What this model shows is that the stability 
of the supercooled liquid should not be taken for granted. 

Finally, even if stability were to be proved not to be a concern in real systems, this remains a cautionary note 
when formulating non-disordered glass models, especially on the lattice, which is so attractive. One can also wonder 
if these considerations might be applicable to disordered systems, or systems with complex (highly degenerate) low- 
temperature phases, were the supercooled liquid could be unstable with respect to a crystal in disguise. 
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